submitted to Journal of Spintronics and Magnetic Nanomaterials 
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Preliminary first-principles calculations on the magnetic behavior of ultra-thin epitaxial multilay- 
ers between the BiFeOs magnetoelectric compound and various types of spacers are presented. As 
spacer we have considered i) InP semiconductor, ii) Fe which is a ferromagnet, and iii) metallic V. 
In all cases under study the growth axis of the multilayer was the [001]. Our results indicate that 
the magnetic properties are seriously downgraded for the ultrathin BiFeOs multilayers independent 
of the nature of the spacer and in some cases under study magnetism even vanishes. More extensive 
calculations are needed to establish a more clear view of the physical properties of the interfaces 
involving the BiFeOs compound. The present manuscript completes the study presented in two re- 
cent research articles [K. Koumpouras and I. Galanakis, J. Magn. Magn. Mater. 323, 2328 (2011); 
ibid, J. Spintron. Magn. Nanomater. 1, in press (2012)]. 
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INTRODUCTION 

Among the most interesting new classes of materi- 
als under intense investigation in Spintronics [IrS| are 
the so-called multiferroics which combine several fcrroic 
orders like ferromagnetism, ferroelectricity, ferroelastic- 
ity ctc.0-@l The compounds combining the electric and 
magnetic order Q have several potential applications like 
magnetic-field sensors and electric-write magnetic-read 
random-access memories. (7- 11 Magnetic order and fer- 
roelectricity have different origins [12Hl4 | and thus the 
materials exhibiting the magnetoelectric effect are few 
and the coupling between the magnetic and electric prop- 
erties is weak. An alternative route to achieve a strong 
coupling could be the growth of thin film heterostruc- 
tures and several advances have been made towards the 
magnetic control of ferroelectricity 
trie control of thin film magnetism. 



15 



18 2 



171 and the elec- 



Bismuth ferrite is probably the most studied represen- 
tative of magnetoelectric materials. Bulk BiFeOs crys- 
tallizes in a perovskite-like pseudocubic structure instead 
of a ferrite one and is classified ferroelectric G- 

type antiferromagnet.[2lJ Several first-principles calcu- 
lations have been carried out to study the pro per ties 
of bulk BiFeOs [22j and we refer readers to Ref. [23| for 
an overview of the literature on this compounds. Since 
the single-component crystals like BiFeOs present only 
a weak magnetoelectric effect, an alternative route to 
achieve a more strong effect has been proposed to be the 
growth of heterostructures where epitaxial strain can en- 
hance the phenomenon. [24 - 27 1 Thus first- principles cal- 
culations of such heterostructures involving alternating 
layers of BiFeOs and various spacers can serve as a test- 
ground to study the behavior of electric and magnetic 
properties of films. The study of the latter in the case 
of ultrathin epitaxial films is the aim of the present 
manuscript. 



In a recent publication (Ref. 1281 ) we have pre- 
sented extended first-principles calculations, employ- 
ing the Quantum-ESPRESSO[29j ab-initio electronic 
structure method in conjunction with the Generalized- 
Gradient Approximation (GGA) in the Perdcw-Burke- 
Erzenhof formulation, (3^] on the electronic and magnetic 
properties of BiFeOs alloy as a function of the lattice con- 
stant in the case of the cubic perovskite structure (see 
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figure 1 in Ref. |28|). In Ref. 

study to cover also the case of Mn substation for Fe 
compounds under study in these two references exhib- 
ited significant magnetic properties with spin magnetic 
moments at the Fe and Mn sites of the order of 3 us- 
In the present contribution we expand these two studies 
using the same ab-initio method with grids of the same 
density in the k-spacc to cover also the case of ultrathin 
epitaxial films. Although these films arc ideal cases and 
cannot be, with few exceptions, be realized experimen- 
tally, they can serve as starting point to understand the 
interplay between the interface structure and the mag- 
netic properties of BiFeOs films with other spacers. For 
completeness we took into account three different spac- 
ers to cover different cases of electronic properties at the 
interface: InP semiconductor, metallic V and ferromag- 
netic Fe. In Section II we present the structure of the 
multilayers and the results for the case of the InP spacer, 
in section III the cases of Fe and V metallic spacers and 
finally in Section IV we present the main conclusions of 
our work. 



InP SPACER 
Structure of the multilayer 

Along the [001] direction BiFeOs can be considered to 
consist of alternating layers of BiO and FCO2 as shown 
in Fig. [TJ The in-plane unit cell is a square of the same 
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FIG. 1: Schematic representation of the structure of the BiO 
(a) and Fe02 (b) possible terminations of the BiFeC>3 alloy 
along the [001] direction. With the red lines we denote the 
limits of the two-dimensional unit cell of InP (see Fig. [2]) 
which has a lattice constant: ai n p = y/2 x aBiFe0 3 - 



lattice constant as the cubic unit cell of the bulk BiFcC>3. 
InP crystallizes in the zincblende structure, similar to 
GaAs, as most of the III-V semiconductors, and has an 
experimental lattice constants of 5.87 A. Along the [001] 
direction, the zincblende structure can be viewed as con- 
sisting of alternating pure In and pure P atomic layers as 
shown in Fig. [2] To fully describe the zincblende struc- 
ture we have to consider, except the In and P atoms, also 
the occurrence of vacant sites (voids) ; there are exactly 
two no-equivalent voids within the zincblende unit cell. 
As shown in Fig. [5] along the [001] direction there are 
not two but four alternating non-equivalent atomic lay- 
ers. The difference between the A(B) and C(D) layers 
is that the In(P) atoms have exchanged sites with the 
voids. This exchange of atomic positions is important 
since it lead to four non-equivalent interfaces between 
InP and BiFc03 in our study. Finally we have to men- 
tion that the lattice constant of InP coincides with the 
BiFcC>3 lattice constant multiplied by the square root of 
two: ai n p x aBiFeOa- Thus as shown in Fig. [1] 

the diagonal of the BiFeOs two-dimensional square unit 
cell can be assumed to be the side of the corresponding 
unit cell of the InP (denoted with red color in the figure) 
and thus epitaxial growth between InP and BiFeC>3 can 
be assumed. 

For the InP/BiFe03 multilayers we studied four dif- 
ferent cases with respect to the relative position of 
the atoms at the interface between the two spac- 
ers. Along the growth axis we took into ac- 
count eight atomics layers which are repeated in the 
[001] direction. Initially we studied the structure 
...Fe0 2 //In/P/In/P/BiO/Fe0 2 /BiO/Fe0 2 //In..., and 
thus in the first case we have two non-equivalent inter- 
face in our structure: a P/Bi and an In/Fe ones. A close 
examination of Figs. [1] and [2] reveals that, in the case 
just mentioned :(a) with respect to the P/Bi interface P 
atoms are situated in the diagonal connecting Bi atoms 
(compare (a) in Fig. [I] with the D layer in Fig. [2]), and 
(2) with respect to the In/Fc interface In atoms are situ- 
ated in the midpoint between Fe nearest neighbors (take 
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FIG. 2: Structure of the atomic layers of InP along the [001] 
direction taking into account also the vacant sites. The struc- 
ture is built using four alternating pure atomic layers A, B, 
C and D. The difference between A(B) and C(D) is that the 
In(P) atoms and the voids have exchanged sites. 



into account A layer in Fig. [2] and combine it with the 
(b) in Fig. |TJ. 

In the second case under study In and P atoms have 
exchanged sites with respect to the hrst case under 
study and thus now we have along the [001] direction a 
. . .Fe0 2 / /P/In/P/In/BiO/Fe0 2 /BiO /Fe0 2 / /P. . . struc- 
ture, i.e. the inequivalent interfaces are now In/Bi and 
Fc/P. Notice that as in the Fe/In interface in case 1, also 
in the P/Fe interface, which occurs in case 2, the P atoms 
at the interface are located in the midpoints between the 
Fe atoms (if we examine only the in-plane projection of 
the multilayer) just above the oxygen atoms and this re- 
sults to vanishing magnetism in case 2 under study as 
will be discussed later in the next subsection. 

As we mentioned above the structure of InP along the 
[001] direction should be viewed as consisting of four 
atomic layers (see Fig. [2]) since voids play a crucial role 
in interfaces. If in the two previous cases 1 and 2 we 
exchange layer A(B) with layer C(D) in Fig. [2] we get 
two new cases 3 and 4 with different arrangement of the 
atoms at the interface. In case 3(4), the succession of the 
atomic layers is similar to case 1(2). When we compare 
case 1(2) with case 3(4), we conclude that the P(In)/Bi 
interfaces are similar since the P(In) atoms are now situ- 
ated at the other diagonal connecting the Bi atoms. Crit- 
ical is the other In(P)/Fe0 2 interface since now in case 
3(4) the In(P) atoms are located just above the Fe atoms 
while in case 1(2) they were located above the oxygen 
atoms at the midpoints between neighboring Fe atoms. 
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FIG. 3: Density of states (DOS) projected on the Fe d-orbitals 
for the Fe atoms in the BiFeC>3 spacer (Fe bulk at the middle 
of the spacer and Fe lnt at the FeC>2 interface for all four cases 
under study (see text for explanation). Positive DOS corre- 
sponds to the majority spin (spin-up) electrons and negative 
DOS to the minority spin (spin-down) electrons. The Fermi 
level has been assigned to the zero energy. 



Results and discussion 

In all cases under study our results converged to the 
ferromagnetic solution independently of the initial con- 
ditions and initial arrangement of the spin moments. As 
in-plane lattice constant we have chosen the \f{2) x 4.153 
A = 5.874 A (11.1 a.u.) which is the experimental lattice 
constant of InP and moreover for the lattice constant 
of 4.153 A BiFc0 3 exhibits very pronounced magnetic 
properties. [11] As a result of the epitaxial growth along 
the [001] axis the lattice constant was 14.182 A (26.8 
a.u.). 

In Fig. |3] we present the density of states (DOS) of 
the Fe d-orbitals for all four cases under study. We de- 
note as Fe bulk the Fe atoms within the BiFeOa spacer 
and with Fe m * the Fe atoms located at the interface. 
We have also to note here that for the multilayer struc- 
tures under study, we have a complete lift of the degen- 
eracy of the Fe d-orbitals and we cannot refer anymore 
to the double-degenerated e g and triple degenerated t^g 
orbitals as in bulk B1FCO3. In case 2 where we have a 
P/Bi interface our calculations have converged to a non- 
magnetic solution. The loss of magnetism in this case 
should be attributed to the reduced hybridization be- 
tween the neighboring Fe and P atoms at the interface. 
The latter have as valence 3p electrons which are less ex- 
tended in space with respect to the In 5p valence states 
and do not hybridize with the Fe d-orbitals at the inter- 
face. In case I, the d-band is shifted to higher energy val- 
ues for the spin-down electrons with respect to the spin- 
up electrons. This is more easy to visualize for the Fe m * 



TABLE I: Fe spin magnetic moments in [lb within the BiFeOs 
spacer for the InP/BiFe03 epitaxial ultrathin multilayers for 
all four cases under study (see text for explanation of different 
cases). The BiFeOa spacer contains two inequivalent atoms: 
Fe bulk at the middle of the BiFe0 3 spacer and Fe lnt at the 
interface. 



Case 1 
Case 2 
Case 3 
Case 4 



0.357 


0.247 
0.296 



1.158 


0.746 
0.919 



atoms since their hybridization with the p-orbitals of In 
leads to more pronounced electronic properties (smaller 
bandwidth with respect to the Ye bulk atoms and thus 
more intense picks). The latter is also reflected on the 
spin magnetic moments presented in Table U where the 
spin magnetic moments are considerable larger for Fe m * 
with respect to Ye hwlk . 

In Fig. |3] we also present the DOS of the Fe atoms for 
the cases 3 and 4 where we have changed the positions of 
the In(P) atoms at the interface. More precisely, as we 
discussed in the previous section, the In(P) atoms at the 
interface with FeO are not any more situated above the 
Fe atoms instead of the O atoms in cases 1 and 2. This 
leads to increased hybridization also in case 4, where we 
have a P/Fe interface and now the multilayer converges 
to a magnetic solution contrary to case 2 where we had 
a non-magnetic configuration. If we compare cases 1 and 
3, where the Fe/In contact appears, we can notice that 
in case 3 the splitting of the Fe mt d-bands is smaller due 
to the alteration in the Fe 3d-In 5p orbitals hybridization 
resulting also in smaller spin magnetic moments as shown 
in Table HI But overall the obtained DOS are similar to 
the ones obtained for case 1 and presented also in Fig. |3] 



V AND Fe SPACERS 

We continue our study using two metallic spacers, fer- 
romagnetic Fe and normal metallic V, instead of InP. 
Both Fe and V crystallize in the bcc structure. The 
zincblende structure of InP, if we take into account the 
vacant sites and ignore the different chemical species, is 
in reality a bcc one. Thus the structure of the Fc/BiFc03 
and V/BiFe03 multilayers is similar to the structure of 
the InP/BiFc03 where all sites are occupied exclusively 
by Fe or V atoms. We have studied 5 cases for both mul- 
tilayers where we have just varied the lattice constant. 
Case 1 corresponds to the same lattice constant as for 
InP/BiFe03 studied in the previous section. Cases 2 and 
3 correspond to a uniform compression of the lattice pa- 
rameter in all directions of 10 and 5 % , respectively, 
and cases 4 and 5 of a uniform expansion of the lattice 
parameter of 5 and 10 % , respectively. Our aim is to 




(a) (b) 

FIG. 4: Schematic representations of the V/BiO (a) and 
V/Fe02 interfaces in the case of the V/BiFe03 multilayer. 
For the Fe/BiFe03 multilayer under study Fe atoms simply 
substitute the V atoms since both V and Fe have a bcc struc- 
ture. 



study the behavior of the magnetic properties upon hyr- 
drostatic pressure. In cases 1, 4 and 5 we converged to 
a ferromagnetic solution irrespectively of the initial ar- 
rangement of the spin magnetic moments. In contrast 
in the compressed cases 2 and 3 for Fe/BiFe03 we could 
not get convergence while for the case of the V spacer 
we converged to a non magnetic solution. This behav- 
ior stems from the reduced volume around the Fe atoms 
which leaded to a large compression of the Fe d-orbitals 
and suppression of magnetism as expected by the mag- 
netovolume effect which has been extensively studied for 
transition metal atoms. In the multilayers under study 
we have also two inequivalent interfaces: a V(Fc)/BiO 
and a V(Fe)/FcC>2 contact. In the first interface the four 
V(Fe) atoms are located within the two diagonals con- 
necting the Bi atoms as shown in Fig. [4jt while in the 
second interface they are located just above the Fe atoms 
of FeC>2 as shown in Fig. |4Jd. 

We will start our discussion concentrating on the Fe 
atoms within the BiFeC>3 spacer in the case of the 
Fc/BiFcC>3 multilayer and in Table ITT1 we have gathered 
the Fe spin magnetic moments for both Fe atoms in the 
inside of the BiFeC>3 spacer (Fe bulk ) and at the interface 
(Fe m *). As mentioned above for the compressed cases 2 
and 3 we were not able to converge to a solution. For 
case 1 we converged to a non-magnetic solution and for 
the more expanded lattice parameters (cases 4 and 5) 
we converged to a magnetic solution. The spin magnetic 
moments are larger for the more expanded case 5. This 
behavior is expected by the well studied magnctovolume 
effect. For the late transition metal atoms like Fe, the 
spin-splitting of the d-states increases with the atomic 
volume. Thus as we expand the lattice the tendency 
to magnetism increases leading to larger spin moments 
while compression of the lattice eventually leads to loss 
of magnetism. A similar situation occurs also for the 
case of V/BiFcC>3 but now the tendency to magnetism is 
stronger for the same lattice parameter. For cases 2 and 
3 we converged to a non-magnetic solution while we got 
a ferromagnetic configuration for case 1 contrary to the 
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FIG. 5: DOS projected on the Fe d-orbitals for the Fe atoms 
within the BiFe0 3 spacer for both V/BiFe0 3 and Fe/BiFe0 3 
multilayers. Results are for the two expanded lattice param- 
eters denoted as cases 4 and 5 (see text for details). Details 
as in Fig. [3] 



TABLE II: Fe spin magnetic moments in /is within the 
BiFe0 3 spacer for the V(Fe)/BiFe0 3 epitaxial ultrathin mul- 
tilayers for all cases under study. Case 1 corresponds to a 
BiFe0 3 in-plane lattice constant of 4.153 A, cases 2 and 3 to 
uniform compression by 10 and 5 % respectively, and cases 
4 and 5 uniform expansion by 5 and 10 % respectively. Zero 
values means that we have converged to a non-magnetic so- 
lution and "— " that we were not able to achieve convergence. 
Details as in Table U 

V/BiFeC-s VFe/BiFe0 3 
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Case 4 


0.269 


1.252 


0.224 


0.620 


Case 5 


0.325 


1.301 


0.339 


0.815 



Fc/BiFc03 multilayer. Fe spin magnetic moments are 
larger for the V spacer but the largest calculated value, 
which we got for Fe mt in case 5 as shown in Table lU 
is 1.3 hb almost half the value in pure bulk Fe. Al- 
though magnetism is present for the cases with the more 
expanded lattice constant the magnetic properties are 
seriously downgraded with respect to pure BiFeOs bulk 
crystals. 

In the case 1 of Fc/BiFc03, presented in the lower 
panel of Fig. [7J both Fe atoms within the BiFeOs spacer 
are non magnetic and the DOS is the same for both spin 
directions, while in the case of V/BiFc03 a small split- 
ting of the d-bands appears in accordance to the spin 
magnetic moments presented in Table [TTJ This splitting 
increases as we move to case 4 and 5 (presented in Fig. 
[5]) following the increase of the spin magnetic moments. 
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FIG. 6: Same as Fig. U but now we have exchanged the V 
atomic layers. In the left (a) panel the V/BiO interface and 
in the right panel (b) the V/Fe02 interface. Notice that with 
respect to Fig. [3] the V atoms are located now exactly at the 
top of the O and Bi atoms. We do not present the V layers 
and the BiO(Fe02) layers in the same square as in Fig. [4] to 
make the atomic positions more clear. 



We can make two remarks with respect to the presented 
DOS. First, for the Fe atoms at the interface (Fe mt ) the 
weight of the d-states is shifted close to the Fermi level 
with respect to the Fe bulk atoms as a consequence of the 
increased hybridization at the interface. Second, if we 
compare the behavior of the same Fe atom for the same 
lattice constant in the two multilayers under study, the 
exchange splitting of the d-bands is larger in the case of 
the V spacer (spin-up states are deeper in energy and 
spin-down states higher in energy) in accordance with 
the larger spin magnetic moments in this case. 
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FIG. 7: DOS projected on the Fe d-orbitals for the Fe atoms 
within the BiFe0 3 spacer for both V/BiFe0 3 and Fe/BiFe0 3 
multilayers. Results are for the lattice constant of 4.153 A. 
In the bottom panel we present the results for the interface 
structure of Fig. [4] and in the upper panel for the interface 
structure in Fig. [6] Details as in Fig. [3] 



is 0.26 fiB, while the Fe bulk moment is 0.11 [is for the 
case of the V spacer and only 0.002 mwg for the case of 
the Fe spacer. Our discussion on the spin magnetic mo- 
ments is reflected also on the DOS presented in the upper 
panel of Fig. [7] where the imbalance between the spin- up 
and spin-down states is very small. The only noticeable 
effect is the shift of the Fe mt d-states between the two 
multilayers; for the V spacer the Fe mt states are more 



concentrated around the Fermi level. The Fe 



bulk 



atoms 



are shielded from the interface due to the surrounding Bi 
and O atoms and the DOS is more similar for both type 
of V and Fe spacers as reflected also on the spin magnetic 
moments. 



Effect of the change of V(Fe) positions at the 
interface 



Behavior of V and Fe atoms within the V(Fe)spacer 



In this subsection we present results again for the 
V/BiFe03 and Fe/BiFe03 multilayers for the lattice con- 
stant of case 5 for which magnetism is more pronounced. 
And thus the in-plane lattice constant is 6.46 A and the 
out of plane 15.61 A . The difference with the multilayers 
studied just above is that we have changed the position 
of the V(Fe) atoms at the interface. As shown in Fig. 
|6]the V(Fe) atom at the interface with FeO are now lo- 
cated just above the oxygen atoms of the BiFeOs spacer 
and not the Fe ones. The change of the local environ- 
ment leads to reduced hybridization of d-orbitals of the 
transition-metal atoms at the interface due to the larger 
Fe-V(Fe) distance and thus to smaller induced spin mag- 
netic moments in the BiFeOs layer. For both V/ and 
Fc/BiFc03 multilayers the Fe mt spin magnetic moment 



In this section we will shortly refer to the magnetic 
properties of the V and Fe atoms within the transition 
metal spacers. There are four transition metal atoms 
within each atomic layer and we have 4 atomic layers 
in our spacer which we count starting from the interface 
with the Fe02 (we denote it as layer 1) and end with the 
atomic layer at the interface with BiO. We will start our 
discussion for the V spacer in the case of the V/BiFc03 
multilayer. V atoms present vanishing spin magnetic mo- 
ments as shown in Table IIIII in all cases under study and 
only in the case 5 (largest lattice constant under study), 
do the V atoms in the layer close to the Fe02 layer show 
a small spin magnetic moment of about 0.13 /is. 

Similar is the situation for the Fe spacer where for 
the two smaller lattice constants wc could not even con- 
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TABLE III: V(Fe) spin magnetic moments in \ib within the 
V(Fe) spacer for the V(Fe)/BiFe03 epitaxial ultrathin multi- 
layers. We do not present results for the compressed cases 2 
and 3 since the V/BiFe03 was found non-magnetic while for 
Fe/BiFe03 we could not get convergence. We have 4 atomic 
layers of V(Fe) spacer and each one contains four V(Fe) atoms 
as shown in Fig. [3] We count the layers starting from the one 
at the V(Fe)/Fe02 interface and finishing at the V(Fe)/BiO 
interface. With "case 5-II" we denote the structure presented 
in Fig. [6] Notice that within each layer there are two inequiv- 
alent with respect to their magnetic properties Fe(V) atoms, 
e.g at the V/FeCh interface presented in Fig. [3]there are the 
V atoms at the diagonals located just above at the Fe atoms 
and the V atoms at the corners and the middle of the square 
unit cell. Near the BiO interface all V(Fe) atoms as shown 
in Fig. [4] have the same nearest neighbors environment and 
for all atoms in layer 4 close to the V(Fe)/BiO interface we 
got the same value of the spin magnetic moment for all V(Fe) 
atoms within the layer. 

V/BiFe0 3 





case 1 


case 4 


case 5 


case 5-II 


Layer 1 


0.016 


-0.010 


-0.003 


-0.003 




0.062 


0.130 


0.095 


0.039 


Layer 2 


-0.005 


0.015 


0.016 


0.007 


Layer 3 


-0.015 


-0.003 


-0.004 


0.015 




0.002 


-0.006 


-0.009 


-0.002 


Layer 4 


-0.005 


0.005 


0.010 


0.003 






Fe/BiFe0 3 






case 1 


case 4 


case 5 


case 5-II 


Layer 1 





-0.026 


-0.129 


-0.110 






-0.040 


-0.174 


0.479 


Layer 2 





0.017 


0.209 


-0.023 


Layer 3 





-0.005 


-0.021 


-0.022 






-0.012 


-0.202 


0.058 


Layer 4 





0.015 


0.116 


-0.010 



verge our calculations. As shown in Table ITTTI only for the 
largest value of the lattice constant we got significant val- 
ues of the Fe spin magnetic moments, which even for this 
case, are considerable smaller than the spin magnetic mo- 
ments of Fe atoms of the BiFeC>3 spacer (~0.2 \ib for Fe 
spacer compared to ~ 0.8 [is for BiFeOa spacer) and are 
about one order of magnitude smaller than in bulk Fe. 

We have also included in Table Hill the results for the 
interface structure of Fig. [6l denoted as " case 5-II" where 
we have changed the positions of the V(Fe) atoms at the 
interface. Overall also in this case the magnetic proper- 
ties of the spacer are not significant. Only in the case 
of the Fe/BiFe03 multilayer half the Fe atoms at layer 
1 (located at the interface with Fe02) present a signifi- 
cant spin magnetic moment of about 0.5 fiB which is still 
much lower than the bulk Fe value. 



CONCLUSIONS 

We expand our study on the magnetic proper- 
ties BiFcOa presented in Refs. |28| and [3l| to 



the case of ultrathin epitaxial multilayers using the 
Quantum-ESPRESSO first-principles electronic struc- 
ture method, [i^] We have studied several cases of these 
ultrathin epitaxial BiFcOs multilayers using different 
types of spacers covering a wide range of electronic ma- 
terials: InP semiconductor, ferromagnetic Fe and metal- 
lic V. Irrespectively of the spacer the magnetic proper- 
ties of the multilayers were found to seriously degrade 
due to the low-symmetry of the film and in some cases 
magnetism vanished completely. Thus ultrathin multi- 
layers of BiFcOs are not suitable for spintronic applica- 
tions. More extended calculations covering more thick 
multilayers and including quantum molecular dynamics 
are needed in the future to determine in more detail the 
properties of the interfaces of B1FCO3 with other alloys. 

Financial support from the University of Patras (K. 
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